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We present measurements of the inclusive production cross sections of the T(1S) bottomonium 
state in pp collisions at y/s = 1.96 TeV. Using the T(IS') — > decay mode for a data sample 

of 159 ± 10 pb~* collected by the D0 detector at the Fermilab Tevatron collider, we determine the 
differential cross sections as a function of the T(1S) transverse momentum for three ranges of the 
T(15) rapidity: < |y T | < 0.6,0.6 < |y T | < 1.2, and 1.2 < \y r \ < 1.8. 

PACS numbers: 13.85.Qk 



Quarkonium production in hadron-hadron collisions 
provides insight into the nature of strong interactions. It 
is a window on the boundary region between perturba- 
tive and non-perturbative QCD. Recent advances in the 
understanding of quarkonium production have been stim- 



ulated by the unexpectedly large cross sections for direct 
Jfvb and vb(2S) production at large transverse momen- 
tum (j>t) measured at the Fermilab Tevatron collider Q]. 

Bottomonium states are produced either promptly or 
indirectly as a result of the decay of a higher mass state 
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, e.g. in a radiative decay such as \b —*■ T(lS)-f. The 
only detailed studies of T(nS) production at the Teva- 
tron have been done by the CDF Collaboration 0, Q in 
the limited T rapidity range of \y r \ < 0.4 at y/s = 1.8 
TeV, where y = | In E ^Z" , E is the T energy, and p z is 
the T momentum parallel to the beam direction. 

Three types of models have been used to de- 
scribe prompt quarkonium formation: the color-singlet 
model |4j , the color-evaporation model (and a follow- 
up soft color interaction model 0), and the color-octet 
model ■ These models of quarkonium formation lead to 
different expectations for the production rates and polar- 
ization of the quarkonium states, yet many of the model 
parameters have to be extracted directly from the data. 
A recent paper 0, Q successfully reproduces the shape of 
the pt distribution of T states produced at Tevatron en- 
ergies by combining separate perturbative approaches for 
the low- and high-p^ regions. The absolute cross section 
is not predicted by these calculations, which are similar 
to the color-evaporation model. 

In this Letter we concentrate on the production of the 
T(1S) state. A precise measurement of the differential 
T(1S) cross section, using the wide rapidity range acces- 
sible by the D0 detector, will provide valuable input to 
the various quarkonium production models. By recon- 
structing the T(15) through its decay T(15) — * [i + \jT , 
we determine production cross sections of the T(IS') as 
a function of its transverse momentum, in three rapid- 
ity ranges: < \y T \ < 0.6, 0.6 < \y r \ < 1.2, and 
1.2<|y T |<1.8. 

The D0 detector is described in detail elsewhere [loj . 
Here, we briefly describe only the detector components 
most relevant to this analysis. The D0 tracking system 
consists of a high-resolution silicon microstrip tracker 
(SMT) surrounded by a central scintillating-fiber tracker 
(CFT) inside a 2 T magnetic field provided by a super- 
conducting solenoid. The tracking volume extends to a 
radius of approximately 52 cm. Closest to the interaction 
region is the SMT with a typical strip pitch of 50-80 /im. 
It has a barrel-disk hybrid structure and provides track- 
ing and vertexing coverage in the pseudorapidity range 
\r]\ < 3.0, where r\ = — ln[tan(#/2)] and 9 is the polar 
angle. The CFT consists of eight concentric cylinders of 
pairs of scintillating-fiber doublets. On each cylinder, the 
inner doublet runs parallel to the beam axis and the outer 
doublet is mounted at a stereo angle of ±3°, alternating 
with each cylinder. Located outside the superconducting 
coil is the uranium-liquid-argon calorimeter. Beyond the 
calorimeter, the muon system consists of three layers of 
drift tubes, 10 cm wide in the central region (|?7| < 1) and 
1 cm in the forward region (1 < \r/\ < 2). Interspersed 
between the drift tubes are scintillating counters. Lo- 
cated between the innermost and the middle layers of 
drift tubes are 1.8 T iron toroid magnets. D0 uses up to 
three levels of triggers to reduce the initial event rate of 
1.7 MHz to an output rate of approximately 50 Hz. 




Dimuon mass (GeV/c 2 ) 



FIG. 1: Example of fits to the dimuon spectra in different 
bins of rapidity in the pr bin of 4 GeV/c < < 6 GeV/c: 
(a) \y T \ < 0.6, (b) 1.2 < |y T | < 1.8. The heavy line shows 
the combined fit for signal and background. Also shown are 
the individual contributions from the three T states and the 
background separately. 



The data were collected between June 2002 and 
September 2003 and correspond to an integrated lumi- 
nosity of 159 ± 10 pb _1 for the chosen two triggers. 
These triggers are scintillator-based dimuon triggers at 
the first trigger level and require the confirmation of one 
or both muons at the second trigger level. The first trig- 
ger level is almost fully efficient for muons with a trans- 
verse momentum above 5 GeV/c. For events passing 
our analysis criteria, the second level trigger requirement 
kept more than 97% of events which satisfied the first 
level requirement. 

The analysis requires two oppositely charged muons 
with pf^ > 3 GeV/c and < 2.2. Only muons that are 
matched to a track found by the central tracking system 
and which have hits inside and outside the toroidal mag- 
nets are used. The track associated with a muon must 
have at least one hit in the SMT. We reject cosmic ray 
muons based on timing information from the muon sys- 
tem scintillators. Compared to muons from the dominant 
bb background, muons from T(nS) decays are expected 
to be relatively isolated, and therefore we require at least 
one of the muons to satisfy the following criterion: cither 
the sum of the transverse momenta of charged tracks in 
a cone of radius 0.5 (in rj-(f) space) around the muon is 
less than 1 GeV or the sum of the calorimeter transverse 
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energies in an annular cone of radii 0.1 and 0.5 around 
the muon is less than 1 GeV. This isolation requirement 
reduces the background by 35% and the signal by less 
than 6%. 

Two typical examples of dimuon mass distributions in 
different rapidity bins are shown in Fig. ^ In each plot a 
strong T(15) signal can be seen, accompanied by a shoul- 
der attributed to unresolved signals due to T(2S) and 
T(35) production. The mass distributions are fit start- 
ing from 7.0, 7.5 or 7.8 GeV/c 2 , depending onpj and y r , 
to 13.0 GeV/c 2 using separate mass resolution functions 
for each of the T(nS) states and a third-order polyno- 
mial for the background. The mass resolution function is 
approximated by a sum of two Gaussians with the rela- 
tive contribution and width of the second Gaussian fixed 
with respect to the first Gaussian. The values of this 

I 



where £ is the integrated luminosity for the data sample 
used, N(T(1S)) is the number of observed T(15), and 
the £ and k represent the various efficiency, acceptance 
and correction factors. The T(IS') acceptance and recon- 
struction efficiency e acc represents the fraction of gener- 
ated T(1S) events that are successfully reconstructed in 
the D0 detector, not taking into account any loss in ef- 
ficiency due to triggering. Its value is based on a Monte 
Carlo analysis. The dimuon trigger efficiency Etrig for 
reconstructed T(IS') events that satisfy our analysis cri- 
teria is estimated using a trigger simulation and verified 
directly with the data using other triggers. The remain- 
ing factors in Eq. ^ account for the differences between 
the data and Monte Carlo and are referred to as correc- 
tions, rather than efficiencies. The correction k qua i takes 
into account differences in the track quality requirements, 
i.e. the isolation and SMT hit requirements and cosmic 
ray rejection. It is consistent with being independent of 
Pt and its value varies between 0.85 and 0.93 with in- 
creasing rapidity. The central tracking correction ktrk 
takes into account both differences in the tracking and 
the track-to-muon matching efficiency. It is derived from 
the 3 1 ip data sample and Monte Carlo simulation and is 
very close to unity except for the forward rapidity region 
where k tr k ~ 0.95. The correction factor k^imu accounts 
for the differences in the local (i.e. muon system only) 
muon reconstruction, taking into account trigger effects. 
It was determined using J/ip candidates collected with 
single muon triggers. It does not show a significant px 
dependence, but it changes with the muon rapidity. 

In TablelJwe summarize the values of efficiencies found 



contribution were determined from Monte Carlo studies 
and J/tp signal fits to data. The mass of the T(1S) is a 
free parameter of the fit and the remaining two masses 
are shifted by the m(Y(nS)) — m(T(15)) differences of 
563 MeV/c 2 (T(25)) and 895 MeV/c 2 (T(35)), taken 
from Ref. 11]. In addition, only the width of the T(LS') 
state is allowed to vary. The widths of the other states 
are assumed to scale with the mass of the resonance. 
Normalizations of functions representing each resonance 
are free parameters of the fit. The Monte Carlo samples 
used in this study were generated with PYTHIA v6.202 
|l2| . The muon kinematic distributions from data and 
Monte Carlo agree within a given p^, and y bin. 

The cross section for a given kinematic range, multi- 
plied by the branching fraction T(IS') — > fi + n~ , is given 
by: 



(1) 

I 

in different rapidity regions. The measured cross sections 
are collected in Table [D] These cross sections are nor- 
malized per unit of rapidity. 



TABLE I: Efficiencies used in the cross section calculations. 



\v 


Sac 




Ztrig 


k qua l 


ktrk 


kdimu 


0.0 - 0.6 


0.15 - 


0.26 


0.70 


0.85 


0.99 


0.85 


0.6 - 1.2 


0.19 - 


0.28 


0.73 


0.85 


0.99 


0.88 


1.2 - 1.8 


0.20 - 


0.27 


0.82 


0.93 


0.95 


0.95 



Differential cross sections, normalized to unity, are 
summarized in Table ITTT1 Figure [3 shows these cross sec- 
tions compared to theoretical predictions from Ref. |IJ. 
There is little variation in the shape of the px distribu- 
tions with rapidity. This is further illustrated in Fig. [3| 
which shows the ratio of the differential cross sections of 
a (1.2 < \y T \ < 1.8) to a(\y T \ < 0.6). In Fig. Hwe show 
a comparison with results from CDF 

The overall systematic uncertainties, excluding lumi- 
nosity, are app roximately 10%. The uncertainty on the 
luminosity [l3j is 6.5%. The main systematic errors are 
due to the fitting procedure and the determination of 
kdimu- The statistical uncertainty of the fitted number 
of events in a given kinematic bin and the uncertainty 
from varying the contribution of the second Gaussian 
are added in quadrature to give the uncertainties labeled 
'stat' in Table |n] The net effect is an increase in the 
overall fit uncertainty by less than 40% of its statistical 
uncertainty alone. An additional uncertainty in the fit- 
ting procedure due to varying the fitting range and the 



^ t ^W(i*)-mV-)- N{T(1S)) 



dpr ■ dy C- Ap T ■ Ay ■ e acc ■ s trig ■ k qual ■ k trk ■ k dlmu 
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TABLE II: Fitted number of events and d<r(T(lS))/dy x B(T(15) 



per unit of rapidity. 



\y 



Number of T(lg) 



d<x(T(lS))/dy (pb) 



0.0 - 0.6 
0.6 - 1.2 
1.2 - 1.8 
0.0 - 1.8 



12,951 ± 336 
16,682 ± 438 
17,884 ± 566 
46,625 ± 939 



732 ± 19 (stat) 
762 ± 20 (stat) 
600 ± 19 (stat) 
695 ± 14 (stat) 



± 73 (syst) ± 48 (lum) 
± 76 (syst) ± 50 (lum) 
± 56 (syst) ± 39 (lum) 
± 68 (syst) ± 45 (lum) 



TABLE III: Normalized differential cross sections for Y(IS') in different rapidity regions. Quoted uncertainties include statistical 
uncertainties added in quadrature to systematic uncertainties due to the assumed shape of the mass resolution function (cf. 
'stat' uncertainties in Table ITT|l . The remaining systematic uncertainties are pr independent and quoted in Table ITT1 



i (GeV/c) 


0.0 < \y r \ < 0.6 


0.6 < \y T \ < 1.2 


1.2<|y T |<1.8 


0.0 < \y r \ < 1.8 


- 1 


0.051 ± 0.005 


0.061 ± 0.006 


0.050 ± 0.005 


0.056 ± 0.004 


1-2 


0.138 ± 0.010 


0.137 ± 0.010 


0.136 ± 0.011 


0.136 ± 0.008 


2-3 


0.152 ± 0.010 


0.153 ± 0.010 


0.175 ± 0.015 


0.160 ± 0.009 


3-4 


0.149 ± 0.011 


0.175 ± 0.012 


0.160 ± 0.014 


0.159 ± 0.009 


4-6 


0.112 ± 0.006 


0.110 ± 0.007 


0.115 ± 0.008 


0.113 ± 0.005 


6-8 


0.067 ± 0.005 


0.061 ± 0.004 


0.056 ± 0.005 


0.062 ± 0.003 


8-10 


0.034 ± 0.003 


0.034 ± 0.003 


0.034 ± 0.003 


0.035 ± 0.002 


10 - 15 


0.014 ± 0.001 


0.011 ± 0.001 


0.011 ± 0.001 


0.012 ± 0.001 


15 - 20 


0.0032 ± 0.0005 


0.0019 ± 0.0003 


0.0019 ± 0.0004 


0.0023 ± 0.0002 



background parameterization is at the 4% level. The sys- 
tematic uncertainty for kdi mu is 8.7%, 8.2% and 7.2% for 
the three rapidity bins. These were derived from uncer- 
tainties for the Monte Carlo - data differences for indi- 
vidual muons, determined as a function of muon rapid- 
ity and transverse momentum. The other uncertainties 
considered include momentum resolution, uncertainties 
introduced by the track quality and track matching re- 
quirements, variations in the input Monte Carlo distri- 
butions, and changes in detector performance over time. 
All these systematic uncertainties contribute less than 
2% each. 

The current analysis assumes that the T(1«S) is pro- 
duced unpolarized, in agreement with the CDF measure- 
ment Q of the polarization parameter a = —0.12 ± 0.22 
for 8 < p^ < 20 GeV/c. Although we do not include a 
contribution to the systematic uncertainty due to this as- 
sumption, we estimate the sensitivity of our results to the 
Y(IS') polarization by varying a within ±0.15 (±0.30). 
This changes our results by less than 4% (15%) in all pt 
bins. 

In conclusion, we present a measurement of the inclu- 
sive production cross section of the T(1S) bottomonium 
state using the T(IS') — * decay mode. The mea- 

sured cross section x B(T(1S) — ► for the |y T | < 

0.6 region is 732 ± 19 (stat) ± 73 (syst) ± 48 (lum) pb. 
Taking into account a predicted increase in the cross sec- 
tion when the pp center-of-mass energy increases from 1.8 
TcV to 1.96 TeV [l2J, our result is compatible with the 
CDF result 0] of 680 ± 15 (stat) ± 18 (syst) ± 26 (lum) 
pb for =1.8 TeV. We measure the ratios of the cross 
sections for the 0.6 < \y r \ < 1.2 and 1.2 < \y r \ < 1.8 



ranges to that for the \y | < 0.6 range to be 1.04 ± 
0.14 and 0.80 ± 0.11, compared with predictions from 
Monte Carlo [13 of 0.94 and 0.83. Between the rapid- 
ity regions, there is little variation in the shapes of the 
differential cross sections, and their shapes agree reason- 
ably well with theoretical predictions jpj. The shape of 
the combined differential cross section for |y T | < 1.8 is 
consistent with the CDF measurement in the limited ra- 
pidity range of |y T | < 0.4 The results presented in 
this Letter will allow a more precise determination of pa- 
rameters of the various bottomonium production models. 
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FIG. 2: Normalized differential cross sections for T(1S) pro- 
duction compared with theory predictions 0, Oj . The errors 
shown correspond to the errors in Table ITTT1 



FIG. 4: Normalized differential cross sections for T(IS') pro- 
duction at y^s = 1.96 TeV compared with published CDF 
results 3] at yfs = 1.8 TeV. The errors shown are statistical 
only. 
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FIG. 3: The ratio of differential cross sections (squares) for 
cr(1.2 < \y r \ < 1.8) to <d]v T \ < 0.6). The solid line is the 
Monte Carlo prediction (1211 normalized to the measured ra- 
tio of the pT-integrated cross section. Uncertainties of the 
relative normalization are indicated by the dashed lines. 
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